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Abstract. This study presents a complete description of the satellite tracking system 
used for Aryabhata. The total system consists of (i) an interferometry system to mea¬ 
sure the direction cosine angles of the satellite with respect to the ground station, 
(ii) a tone-ranging system to measure the range of the satellite and (iii) a Doppler 
system to measure the range rate of the satellite. The system accuracies obtained 
were ±10' of an arc in the direction cosine angles, ±T0 km in the range and ±6 - 6 
m/s in the range rate. 
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Doppler system. 


1. Introduction 

The tracking system is an essential part of a ground station for obtaining the orbital 
elements of any satellite. Since orbit degradation occurs under perturbing forces 
such as drag due to atmospheric density, earth oblateness etc., it is necessary to 
continuously track the satellite in order to update the orbital parameters. 

Theoretically, six independent observations of a particular tracking parameter of 
a satellite are sufficient to determine its orbit completely. However, due to inaccura¬ 
cies of measurement involved, more than one tracking parameter is usually measured 
and weighting for the data of a particular parameter is given depending upon the 
accuracy of measurement, for orbital computation purposes. 

For Aryabhata, because of spacecraft constraints, no separate onboard tracking 
package was installed. However, with the available onboard communication 
packages (telecommand receiver and telemetry transmitter), three tracking systems 
were developed and are described below. 

(i) An interferometry system, to provide angular information of the satellite. It 
uses the carrier component of the telemetry transmission as the beacon for 
tracking purposes. 

(ii) A tone-ranging system, to provide range. It uses the onboard communication 
packages as transponder for the measurement of range. 

(iii) A Doppler system, to provide range rate. It uses the onboard telemetry 
transmitter carrier frequency as the beacon frequency for tracking puiposes. 

The three tracking systems were installed at SHAR and are being used to obtain 
the tracking parameters of Aryabhata having a near circular orbit of 600 km altitude 
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at an inclination of 51° to the equatorial plane. The station at SHAR provides 
tracking data for about 10 min for the satellite pass with a maximum elevation of 90°. 
The data are recorded using a digital printer. 

The tracking data are provided to the Post-Launch Operation Group for analysis. 
The analysis consists of smoothing the data by a curve fitting procedure and then 
using the smoothed data to up-date the orbital parameters of the satellite. The 
following sections deal with the requirements, principles, design considerations and 
performance of the tracking systems. 


2. Choice of the systems 

The systems conceived for Aryabhata are relatively simple compared to available 
systems due to the limitations of spacecraft configuration and availability of resources 
and time. The three tracking systems—interferometry, tone-ranging and Doppler 
tracking, the last being a supplement to the first two—were decided upon to meet the 
following requirements of accuracy: 

direction cosine angle : ±6' of an arc; 

range : ±1 km; 

range rate : ±6-6 m/s. 

Sophisticated interferometry systems capable of giving an accuracy of better than 
±20" of arc in the direction cosine angles have been developed over the years by 
many nations. This involves huge antenna structures providing a directive gain, 
automatic electrical length measurement and control systems and an elaborate 
calibration procedure (Nollet et al 1974). 

For Aryabhata, as the requirement of accuracy was only ±6' of arc, it was 
decided to go in for a simple interferometry system consisting of six omnidirectional 
(turnstile) antennas with highly phase-stable foamflex cables to feed four pairs of 
receivers located in the operator’s room. 

The analysis indicates that with this simple system it is possible to meet the require¬ 
ment of the angular accuracy. 

There are three widely used systems for the measurement of the range of a space¬ 
craft from the gound station: 

(i) a coded ranging system, 

(ii) a multiple continuous-wave-tone ranging system, 

(iii) a hybrid of the above two. 

These systems differ from each other in accuracy, acquisition time, maximum range 
capability, and compatibility with other signals to be used with the ranging signal. 

A coded ranging system cannot be used with the PCM telemetry signal of interest. 
So, the choice is mainly between a multiple CW-tone ranging system and a hybrid 
system. Analysis shows that the time required to acquire the hybrid signal is 
inversely proportional to the signal-to-noise ratio (SNR) at the input, while in the 
case of multiple CW-tone ranging system, it is independent of SNR. Though a 
hybrid system can give better accuracy, it is more complex than a CW-tone ranging 
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system. Considering these, the latter was selected for ranging purposes (Kuwand 
et al 1974; Peterson and Gupta 1965; Baghdady and Kruse 1965). 

A two-way Doppler system with an onboard coherent transponder provides 
better accuracies than a single way Doppler system, as onboard oscillator instability is 
cancelled. However, for Aryabhata, a one-way Doppler system was selected because 
a separate transponder could not be installed onboard due to various spacecraft 
constraints like power, size and weight, and ground station constraints like 
availability of resources, time and complexities involved. 


3. Interferometry tracking system 

3.1. The principle 

This system is used to give the direction cosine of the satellite with respect to north- 
south and east-west directions. In this system there are basically two receiving 
antennas separated by a distance D (figure 1). For a distant target carrying a 
beacon, the direction cosine p with respect to the line connecting the antenna is given 
by 


p = COS a — j, A / D, 


( 1 ) 


where A is wavelength of transmission, a is direction cosine angle, and <f> is phase 
difference of the wave arriving at the two antennae expressed in a fraction of a cycle. 
By partially differentiating (1), we get 


and 


dp_dh djDd$ 

p T A D + ’ 


da = cot a 


dh _ dD Qj>' 
A D ^ _ 


( 2 ) 

(3) 


The sign of quantities in brackets indicates only the direction of errors involved 
and in the worst case they become additive. Hence, for the worst case, 


da = cot a 


[ 


5A QD 8<f> 
T ^ ~D + J 


]' 


(4) 



Figure 1. Diagram showing the basic principle of the interferometry system 
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It can be seen from (2) that the accuracy in the determination of p can be increased 
by increasing D for a given error in the measurement of <f>. But if D is increased 
beyond A, there is an ambiguity in the determination of p, since <j> is cyclic for a path 
difference of more than A. Hence, generally two or three pairs are provided in each 
direction so that coarse, medium and fine accuiacies are obtained and the ambiguity 
in the fine channel can be resolved using medium and coarse channels. 

3.2. System design considerations 

The system conceived for Aryabhata has two pairs of antennae, in each of the two 
orthogonal directions (N-S and E-W), one as a coarse channel and the other as a 
fine channel. The following achievable errors were assumed before determining 
the base-length. 

(i) An error <9</> of ±7° in the measurement of phase difference, arising from 
differential phase stability of the circuit and the phase jitter for a signal level 
of —130 dBm. 

(ii) An error fiA in the assumption of the wavelength of transmission amounting 
to the frequency instability of the onboard transmitter plus Doppler shift in 
frequency (±3-5 kHz). 

(iii) An error dD in the determination of the distance between the antennae 
(±2 mm). 

Substituting in (4) the above errors and the requirement of accuracy in a of ±6' 
of arc for 30° <a<90°, we found that a baseline separation of A for the coarse 
channel and 25 A for the fine channel are required. Moreover, with these errors in the 
coarse channel, it is possible to resolve the ambiguity in the fine channel (appendix 1). 



Figure 2. Geometry of the antenna field for the interferometry system 
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3.3. Antennae 

The antenna field essentially consists of 6 antennae forming a total of 4 pairs of 
antennae (figure 2). One antenna each in N-S and E-W directions is used as the 
reference for both coarse and fine channels. Two dummy antennae are added for 
symmetry. The antennae are of the turnstile type and are designed to have nearly 
omnidirectional characteristics over a range of elevation angles between 30° and 90°. 
The calculations were made for various values of dj A, where 4/2 is the height of the 
antenna over a reflective platform and are given in table 1. From the table it is 


Table 1. Tracking antenna gain over isotropic antenna 


4/A -* 


0-5 

0-7 



0-9 

Angle (deg.) 

P 

E 

P 

E 

P 

E 

from zenith 

dB 

dB 

dB 

dB 

dB 

dB 

0 

+4-14 

+4-15 

+ 3-81 

+ 3-81 

-5-63 

-5-63 

15 

+4-11 

+4-03 

+ 3-91 

+ 3-83 

-3-53 

-3-61 

30 

+3-99 

+ 3-93 

+ 3-95 

+ 3-09 

+0-59 

-0-47 

45 

+3-19 

+0-47 

+4.15 

+ 1-43 

+ 3-61 

+0-89 

60 

+ 1T5 

-3-79 

+ 3-15 

-1-79 

+4-35 

-0-59 

75 

+3-89 

-12-41 

-1-19 

-9-71 

+ 0-95 

-7.57 

p = pattern gain 







E = effective signal level including range variation 






figure 3. Block schematic of the interferometry system for one channel 
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seen that for a df A=0-7, the required gain pattern is achievable ensuring that 
once the satellite is acquired, the signal does not fall for higher elevation angles 
and the signal strength is fairly uniform. 

The block schematic of the entire system is shown in figure 3. The special local 
oscillator (SLO) provides two frequencies separated by 125 Hz and phase-locked to 
each other. The signals from a pair of antennae are added and amplified in a 
common channel, to keep the differential phase shift at a minimum. Using this 
scheme, the phase difference which exists at the inputs of the pair of antennae is 
transferred to 125 Hz signal appearing at the output of the square law detector, 
relative to the reference 125 Hz . In this scheme, the error due to differential phase 
shift between the channels is minimised. 


3.4. Specifications 

Antennae 
Mounting 
Base-length 
Pattern gain 
Polarisation 
Locking sensitivity 
Noise figure 
Predetection BW 

Phase lock loop noise bandwidth 
(at 125 Hz) 

Resolution of phase measurement 


: turnstile-nearly omnidirectional 
: non-steerable 
: A and 25 A at 137-44 MHz 
: 3-8 dB to -1-2 dB; 90° >a>30° 
: circular 
: -130 dBm 
: 3 dB 
: 15 kHz 

: 2-5 Hz 
: 0 01 ° 


4. Tone-ranging system 
4.1. The principle 

In this system, the range of the satellite from the ground station is measured by 
utilising the principle that an electromagnetic wave when propagated over a 
distance experiences a delay. If the two way range is 2 r and A is the wavelength, 
the total phase delay suffered by the wave, expressed as a fraction of a cycle, 
is given by 

<{> + n = 2r/A, (5) 

where n is an integral multiple of 2rr radians of phase shift suffered by the wave. 
From (5) it can be noticed that the maximum unambiguous range obtainable depends 
upon the frequency of the tone used for ranging; the lower the frequency, the higher 
is the unambiguous range. 
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Figure 4. Block schematic of the tone ranging system 


4.2. System design considerations 

For a given inaccuracy in the measurement of <j>, increasing the tone frequency reduces 
the error in r. For an unambiguous range of 4687-5 km and an accuracy of ±1 km, 
four tones, 32 Hz, 160 Hz, 800 Hz, and 4000 Hz were selected. The required accuracy 
is obtained at the highest tone and the ambiguity is resolved successively by the 
lower tones (Habib et al 1964; Shaffer et al 1964). The lower tones, before they are 
transmitted are up-converted to 4032 Hz, 4160 Hz and 4800 Hz and added to 4000 Hz 
to conserve the bandwidth. These tones are transponded by the telecommand receiver 
and telemetry transmitter onboard. The ground transmitting section consists of a 
tone generator, an up-converter of lower tones and an AM transmitter. The ground 
receiving section consists of an input filter, tone filters, phase correcting networks, a 
digital processor, display and an interface to printer. The block schematic of the 
tone ranging system is given in figure 4. 

The received 4 kHz tone is extracted from the composite tones using a phase locked 
loop. The phase correction networks are used for correcting any offset in the phase 
delay. With the help of transmitted and received tones, the digital processor generates 
a pulse whose width is the measure of range. This pulse is fed to a counter as a 
gate to give the range directly. 


4.3. System specifications 

Maximum range 
Tone frequencies 

Phase jitter allowed in lower tones 
Phase jitter in 4 kHz tone 
Accuracy in the range measurement 
Uplink 
Downlink 


4687.5 km 

32 Hz, 160 Hz, 800 Hz and 4000 Hz 
±30° (3a value) 

±10° (3(7 value) 

± 1 -0 km (3(7 value) 

AM 

PM (deviation=l rad.) 
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5. Doppler tracking system 

5.1. The principle 

The range rate of the satellite with respect to the ground station is determined by 
the Doppler technique. The Doppler frequency is given by f i =V fjC if C, 

where V is the radial velocity of the satellite, f„ is the transmitter frequency and C 
is the velocity of light. Thus, to determine the Doppler frequency it is necessary that 
the transmitter frequency is known precisely. 

5.2. System design considerations 

The telemetry carrier working at 13744 MHz is used as the beacon signal. This 
necessitated a highly stable frequency source for telemetry transmission with a long 
term stability of 1 x 10' 6 and a short teim stability of 1 X 10' 8 over the average radio 
visibility duration of 10 min. At this carrier frequency, the Doppler shift ranges 
between ±4 kHz for a circular orbit of 600 km altitude. 

The ground segment (figure 5) consists of an antenna with a preamplifier and a 
double superheterodyne receiver. The second IF is at a frequency of 20 kHz which 
also acts as a bias frequency for the Doppler to swing in either direction by ±4 kHz. 
Since any instability in the local oscillators used in the receiver affects the accuracy 



I F. = Intermediate frequency 

Figure 5. Block schematic of the Doppler tracking system 
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Figure 6. Block diagram for the phase lock loop 

of the system, oven controlled crystal oscillators with a stability of 1 in 10 10 were 
used. 

The heart of the system is a tracking filter (a phase locked loop (PLL) operating 
at 20 kHz) providing a small bandwidth of less than 15 Hz to track the incoming 
signal frequency even under the worst SNR conditions. The SNR at the input of 
the receiver is of the order of 0 dB. The PLL filter keeps searching for the incoming 
frequency in the band 20 kHz±4 kHz and when the signal is present it locks. 
Once the PLL filter gets locked to the signal it keeps tracking the frequency changes 
in the incoming signal. The schematic diagram of the Doppler tracking filter is given 
in figure 6. In a second order loop with F{S)= l-f a/S, and with a linearly varying 
frequency as the input signal (as is the case with Doppler) the phase error function 
is given by 


4>(S) = 


S 2 




+ 


W~w 0 e Q- 

s 2 sy 


( 6 ) 


where /l 2 =signal power, /f—loop gain, R=rate of change of frequency in radians/s 2 , 
.S'^complex frequency, W— lT 0 =difference between the centre frequency of VCO 
and input frequency and 0 o =initial phase of the signal. So the steady state phase 
error is 


lim 

/ -*■ oo 


<M0 = 


R 

aAK' 


( 7 ) 


and the loop is in perfect frequency lock. 

For a second order loop, when the received signal frequency varies linearly with 
time, the loop acquires the signal if the frequency slope of the Doppler shifted signal 
is less than half the value of aAK and it can keep track of the signal till the ratio 
becomes unity (Viterbi 1966). The mean square phase jitter is inversely propor¬ 
tional to the loop SNR and is given by 

W =--- ( 8 ) 

"° 2(SNR) l 
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As the frequency is measured by 6ft), there is a frequency error when measurement 
is made over an integration time T t . The error in /is 


A /= 


1 

4jt7’ < (SNR) £ * 


5.3. Specifications 


Sensitivity 
Antenna 
IF bandwidth 
PLL bandwidth 


-130 dBm 

Same as interferometry 
10 kHz 
15 Hz 


( 9 ) 


6. Performance 

6.1. Interferometry 

The interferometry system became functional two weeks after the launch of 
Aryabhata. The spin of the satellite around 50 rev/min caused amplitude 
fluctuations in signal strength. This caused an increase in phase jitter necessitating 
a smoothing of the data before putting them into actual use. Fixed offsets 
in the data were removed after analysis of the data pertaining to various orbits. 
An accuracy of ±10' of an arc was achieved against the design target of ±6' 
of arc. One representative curve of the observed data is given in figure 7. Apart 
from the errors in the system there are many other sources of error such as 
ionospheric refraction of the order of 0-05° to 0 2° at 45° elevation under average 
ionospheric conditions, antenna misalignment in N-S and E-W directions of the 
order of 20" of arc, and ground levelling and ground reflections of the signal. 
However, corrections based on the predictions of the index of refraction can be 
applied. In the initial phase, the entire data were used to generate coarse orbital 
elements but once they became better defined, the data pertaining only to near 
zenith were used to refine the orbital elements. In this zone, the errors due to the 
ionosphere can be minimised to less than 20” of an arc. During installation, 
the distance between the centres of the antennae of the N-S and E-W lines was deter¬ 
mined using tellurometers and theodolites. In practice, the problems of installing 
the system are greater than those of developing it. 

Some improvements in the installed system have been contemplated to achieve 
better accuracy. These include provision of better antenna gain, reduction of the 
differential phase shifts and a better calibration procedure using an aircraft with a 
flash and a beacon. 

6.2. Tone-ranging system 

The tone-ranging system was calibrated with the satellite for zero distance. The 
system was put into operation 15 days after launching. The tone-ranging data are 
straightforward and can be used as such, as the print-out gives the range in kilometres. 




range (km) 
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Figure 7. Comparison between the observed and predicted phase angles for 
the interferometry system 
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Figure 8. Comparison between the observed and predicted range for the 
tone-ranging system 
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Figure 9. Basic geometry showing different para¬ 
meters for a satellite in a circular orbit around 
the stationary earth 


To remove any dispersion in the data, the data were normally smoothed. A system 
accuracy of ±1-0 km was achieved. The data for one representative orbit, namely 
orbit No. 539 from SHAR, are shown in figure 8. The errors in the data consist of the 
system error resulting from the phase jitter of the fine tone, ionospheric refraction 
error, and drift in calibration due to the temperature fluctuations onboard. 

6.3. Doppler system 

The Doppler system was put into operation from the first visible orbit over SHAR. 
The system errors were estimated before it came into operation. The errors in the 
Doppler data consist of local oscillator stability of 1 in 10~ 10 and the frequency error 
arising from phase jitter in PLL. In addition to this, short term stability of the on¬ 
board transmitter (+1 Hz) and ionospheric effects (upto ±5 m/s) under average 
ionospheric conditions at an elevation angle of 45° in VHF band contribute to the 
total error in the Doppler data. 

For a circular orbit at an altitude h and for stationary earth (figure 9) the range 
rate is given by 


r 


vR . vt 

= — cos a sin- 

r (R + h ) 


( 10 ) 


where v— speed of the satellite, R=radius of the earth, a-— angle subtended at the 
centre of the earth by the station and the satellite at closest approach to the station, 
/•—range of the satellite from the station and t=time of elapse from time of closest 
approach. 

The Doppler curve is symmetrical about the time of closest approach of the satellite. 
However, even for a spinning earth, the satellites with a circular orbit give a sym¬ 
metrical curve to a first approximation. Once the point of inflection and its maxi¬ 
mum slope are found by using mathematical techniques, it is possible to determine 
the exact onboard transmitter frequency and hence Doppler frequency versus time. 
The mathematical technique for a single way Doppler system consists of curve 
fitting and the determination of the time of closest approach, from the available 
data of about 10 min. The typical data, observed on orbit number 781 over 
SHAR, are given in figure 10. 
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4 5 6 7 8 9 10 11 12 13 14 15 

time (min) 

Figure 10. Comparison between the observed and predicted range rate for the 
Dopier tracking system 

Table 2. Tracking systems accuracies 


Parameter 

Systems 

Design goal 
(3 a) 

Achieved accuracy 
(3 a) 

Direction 
cosine angle 

Interferometry 

±6' of arc 

±10' of arc 

Range 

Tone ranging 

±1 km 

±1 km 

Range rate 

Doppler 

±6-6 m/s 

±6-6 m/s 


7. Conclusions 

The satellite tracking systems installed at SHAR, which are the first of their kind 
in India, have provided tracking data which have been used to generate satellite 
ephemeris for Aryabhata. Table 2 indicates the achieved tracking accuracies against 
the design goals. 

It should be emphasised that for improved accuracy in tracking by electronic 
systems, the frequency of operation should be above 1 GHz. In addition a good 
calibration system for interferometry is essential. 
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Appendix 1 

The error in the direction cosine angle is given by 


da = cot a 


&X , 3D d<t> 
T h d ^ T - -’ 


Coarse chaxmel: D/A — 1, 

da = cot a (2-1 x 10-5+90 x 10-5+1942 x 10- 5 ). 
Fine channel: D/A = 25, 

da = cot a (2-1 X 10-5+3-6 X 10-5+77 X lO" 5 ). 


The accuracy of the system is given by the fine channel. With the error in the coarse 
channel, it should be possible to resolve the ambiguity between the two adjacent 
ambiguous direction cosine angles in D— 25A channel. 


Table 3. Error calculation for various angles of arrival a 


a (in degrees) 

da 

Coarse channel 

Fine channel 

89 

± 1-1° 

±2-6° 

70 

± 1-2° 

±3' 

50 

± 1'5° 

± 3-5' 

30 

± 2-2° 

±6' 

10 

± 5-2° 

± 15' 


a (in degrees) 


Interval between adjacent 
ambiguous angles 


10 

± 10° 

30 

±4° 

50 

± 3° 

70 

± 2-4‘ 

90 (zenith) 

±2-3‘ 
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